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Abstract 



^ \ The universal, scaled order parameter profiles P±{z/£,) for critical adsorption 

s, 

, of a fluid or fluid mixture onto a wall or interface, and for the extraordinary 

C , 

Q , transition of the semi-infinite Ising model, are discussed theoretically, where 

CJ 

z is the distance from the interface, ^(T) is the bulk correlation length, and 
^ ' the subscript + (— ) refers to the approach from above (below) T^. Recent 

results to first order in the e = A — d expansion are extrapolated to d = 3 
space dimensions and compared with new Monte Carlo results. In order to 
obtain meaningful extrapolations it is crucial that both the exponential de- 
cay at large C, as well as the known algebraic behavior P±{C,) ~ Q"^!^ at 
small C be correctly reproduced. To this end a recently developed novel RG 
scheme involving a z dependent amplitude renormalization is used. Reason- 
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able agreement of our extrapolations with the Monte Carlo results and some 
experimental results is obtained. 



1. INTRODUCTION 

The phenomenon of critical adsorption of fluids or binary mixed fluids onto walls or 
interfaces has attracted considerable attention, both theoretically and experimentally, in 
recent years 0-^,|l6|. It was originally predicted by Fisher and de Gennes Considering 



a semi-infinite fluid bounded by a planar wall at z = and invoking scaling ideas, these 
authors concluded that the phenomenon has two important signatures: 

(i) Between a microscopic scale a and the macroscopic bulk correlation length ^ there 
exists a regime a ^ 2; ^ ^ of distances z from the interface in which the order parameter 
m{z) = (0(x)) at position x = (xj|, z) decays only algebraically ~ z'^/" . 

(ii) The excess order parameter 

/■oo 

= I [ra{z) — m(oo)] dz , (1) 

i.e., the total amount of adsorbed fluid, behaves as nig ~ t~^'^~^^ as r = (T — T^/Tc 0, 
where is the bulk critical temperature. Since v — f3':^ 0.3 in three dimensions, diverges. 

On a qualitative level both signatures (i) and (ii) have been verified experimentally. 
However, with regard to detailed quantitative investigation of the phenomenon and its quan- 
titative comparison with theoretical predictions, much remains to be done. A detailed theo- 
retical analysis of optical data for critical adsorption has been made by Liu and Fisher [Q, 
who found evidence both for (i) as well as for crossover of m(z, r > 0) from a power-law 
regime at small C, = z/^ to exponential behavior, m[z) ~ e"*', at large C- Although they 
were able to extract also some quantitative information about the nature of the order pa- 
rameter profile from the data, the precision and reliability of their estimates was somewhat 
limited, partly because of the limited accuracy of the available experimental data, but partly 
also because of the lack of quantitatively reliable theoretical results needed as input for the 
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analysis of the experiments. In particular, it became clear that detailed quantitative results 
for the scaling functions P±{C) governing the asymptotic form 

m{z) ^ M^lrf P±{z/^) (2) 

for r ^ 0^ would be highly desirable. Here M_ is a nonuniversal constant, which we fix such 
that P_(oo) = 1. The functions P± also depend on the precise definition of the correlation 
length ^, i.e., on the choice of the nonuniversal coefficients in the asymptotic expression 
^(r — s> 0^) ~ ^o" I''"!"'"- We fix these by taking ^ to be the true correlation length, defined 
through the requirement that the bulk two-point cumulant function decays oc e~'^~^''/^ as 
|x — x'l oo. 

Recently two of us [^,0 succeeded in computing the first two terms in the e = A — d 
expansion of P±{z/C,). Upon extrapolation to = 3 space dimensions explicit results for 
P±{C) were obtained, which had the correct physical features, displaying, in particular, 
a smooth and monotonous crossover from the behavior P±{C) ~ C~^^'^ as C ^ to the 
exponential decay at large (. In the present paper we elaborate on these results, extending 
them and checking them versus Monte Carlo results. We begin with a brief summary of the 
e-expansion results of Ref . , recalling that their naive extrapolation to = 3 would yield 
profiles P±{C) with a totally unacceptable, incorrect short-distance behavior. We show that 
this problem can be overcome in a systematic manner by means of a recently developed 
novel renormalization group (RG) scheme ||^, which ensures the proper, exponentiated form 
of the leading short- distance singularity. We then compare the so-obtained extrapolations 
-P±(C) d = 3) with the results of Monte Carlo calculations. Finally, we present the e expansion 
of a number of universal amplitude ratios, estimate their values for d = 3, and compare them 
with experimental estimates, if possible. 

2. RESULTS OF THE e EXPANSION AND BEYOND 

The analytic results of Ref. were obtained by means of RG improved perturbation 
theory applied to the semi- infinite 0^ model Since this model and its physics is reviewed 

3 



in a separate contribution by Diehl in these proceedings, we will focus on the explanation 
of the results. Critical adsorption (or the so-called normal transition) is described by a 
fixed point P*a^; the extraordinary transition by a fixed point V*^. Although being located 
in different regions of parameter space of the model, these fixed points are equivalent in the 
sense that they yield identical results for the scaling functions P± and other quantities to 
arbitrary order of the e expansion. In Ref. ^ the first two terms in the expansion 

P±{C; e) = P±(C; e = 0) + P±(C; e = 0) + 0{e') (3) 

were obtained. Note that in contrast to our convention here the second-moment definition 
of ^ was used in Ref. 0. We denote this latter correlation length as ^, reserving the symbol 
^ exclusively for the true correlation length. The difference between ^ and ^ is fairly small 
near criticality since the associated amplitudes and differ by less than a few per cent 
T0| . The scaling functions P± and P± corresponding to these two conventions can be easily 



transformed into each other using 

p±(C) = p±(ceo^/t). (4) 

Since ^^/i^ = l + 0{e^) and ^o/io = 1 + (e/12)(ll/2 - vrVS) +0(e2) 0, the e expansions 
of P+ and -P+ agree to order e, but those of P_ and P_ differ already at order e. 

As C ^ 0, the two terms in (H) behave as P±(C; e = 0) ^ c± and 9^ P±(C; e = 0) ^ 
(c-t /6() \n(. If one extrapolated this result naively to = 3 by setting e = 1, one would 
obtain a totally unacceptable short-distance behavior of the form P±(C) ~ [1 + (e/6) In^]/^. 
As discussed in Refs. and 0, the correct limiting form is 

P±{( ^ 0) ^ C-^/^(c± + a± C/-' + a'^ C'^ + 6± + •••) • (5) 

Equation (|^) is consistent with in that the e expansion of the latter agrees with the 
limiting form of (0) for small C,. The coefficients have the expansions 

a+ = V2 (-i + (1 - 6Cs - 6 ln2)) + 0{e^) , (6) 
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a'+ = f + 0(6) 
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+ 0(e) 



^ + 0(6) 



&- = -M + 0(e), 



c+ = V2 \l + ^ {6Ce + 2 \n2 - 13) 



c_ = 2 + I (6 - 16 + vrTs) + 0{e^ 



(7) 
(8) 
(9) 
(10) 

(11) 



where Ce = 0.5772 ... is Euler's constant. Due to the difference between ^ and ^, the 0(e) 
terms of a_ and c_ differ from those of d_ and c_ given in Ref. [^]. 

For large arguments the scahng functions have the expected exponentially decaying forms 

P+(C ^ oo) ^ 2v^ {l + e [i + I In 2 - fyr (l - 2v^) + ^ In f^] + 0(6^)} e''^ (12) 



and 



P_(C ^ oo) ^ 1 + [2 - e I (4 - 1/^3)1 + 0{e^)] e'^ 



(13) 

where it should be noted that would get replaced by e"™"*^* ^ for any choice of the 
correlation length other than the true one. 

In light of the discussion given in Ref. we also determined the next-to-leading terms 
in the asymptotic expansion 



n(C) = e-< + Pi's e-'< + Pi-^ e-'^ + .... 
The value of P^i can be read off of (0). For the remaining coefficients one has 



(14) 



,(oo) 



0{e') 



(15) 



and 



p(oo) 
-M-,3 



1 + e 



+ i ln2 - 27r (1 - 2v^) - 4v^ln + 0(e 



(16) 



The above results provide valuable quantitative information about the nature of the 
scaling functions beyond what is known from scaling considerations. The limiting forms of 
P±{C) at large and small ( can be extrapolated to = 3 in a straightforward manner either 
by simply setting e = 1 in the corresponding coefficients or, preferably, by using improved 
extrapolation techniques (e.g., Pade analyses, into which the exact results for d = 2 have 
been incorporated |]12|). To obtain extrapolations P±{(,d = 3) giving a proper description 
of the ( dependence for all values of C is a more challenging problem. A prime requirement is 
that both the leading asymptotic terms at small and large ( be correctly reproduced. As we 
have seen above, the most naive extrapolation procedure — setting e = 1 in @ — grossly 
fails in this regard. A number of alternative extrapolation schemes have been suggested 
T3|,|6[| , but all of these contain some degree of arbitrariness. Only recently a more systematic 



way of handling this problem was developed [§. 

The crux of this method is a specially adapted, novel RG scheme. Its details are beyond 
the scope of the present paper, so we restrict ourselves to a few remarks. In renormalized 
field theory the amplitude of the order parameter 0(x) usually is reparametrized in a z 

— 1/2 

independent fashion to define the renormalized density (f)^ = 0. Here the amplitude 

renormalization factor Z^{u) depends on m, the dimensionless renormalized coupling con- 
stant, but not on z. The basic new element of the approach of Ref. |^ is that in place of 
a, z dependent amplitude renormalization factor Z{u, zfx) is used, where fi is an arbitrary 
reference momentum. This generalized renormalization factor Z is required to absorb in ad- 
dition to the usual ultraviolet singularities (corresponding to poles in e in the dimensionally 
regularized theory) also the short- distance singularities ~ z~^ln{^z). That is, we require 
that the renormalized profile TnR{z, u, r; /x) = Z^^l"^ m(z) have a finite limit 

\im zmji{z,u,T] fi) . (17) 

z — >0 

At one-loop order a convenient choice is 

Z=l + lu\n{l + j^)+0{e). (18) 
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The advantage of this scheme becomes clear when the resulting RG equations are ex- 
ploited. Due to the z dependence of Z, these are somewhat more complicated than usual, 
albeit solvable by standard methods. The resulting RG-improved perturbation theory yields 
scaling functions -P±(C;e) which exhibit for all e > a smooth crossover from the correct 
short- distance form ~ C,^^^^ to the exponential decay at large C^. Accordingly, meaningful 
extrapolations are obtained upon setting e = 1 (see Figs. 1 and 2). 



3. MONTE CARLO RESULTS 

Monte Carlo calculations have proved to be a very powerful alternative tool for studying 
critical phenomena in confined geometries. To get independent information on the scaling 
functions P± we analyzed data of a Monte Carlo simulation of an L x L x D (L = 128, D = 
80) simple-cubic Ising film with periodic boundary conditions in the L-directions and free 
boundaries in the D-direction. The spins were assumed to interact via ferromagnetic nearest- 
neighbor exchange couplings that take the values Ji and J for all bonds lying entirely within 
the two boundary layers and all other bonds, respectively. All bulk and surface magnetic 
fields were set to zero. To observe the universality class of the extraordinary transition, the 
surface coupling Ji was chosen to be supercritical, i.e., larger than the critical value for the 
occurence of a special transition in the thermodynamic limit L, D ^ oo (cf. Refs. Ill^j7|j9[]). 

To model critical adsorption of fiuids it would seem more realistic to take Ji subcritical 
and to add a surface magnetic field hi > (accounting for the wall-fiuid interaction), i.e., 
to consider the so-called normal surface transition of the model. However, as discussed 
in Ref. [0 and demonstrated exactly in Ref. [1^, the surface critical behavior at critical 



adsorption and the normal transition is representative of the same universality class as the 
extraordinary transition. 

The Monte Carlo data show that the behavior of the order parameter mlz) foTz<^^ is 
well described by m{z) oc (z + Ze)~^^'^, where the exponent is in conformity with the value 
= 0.519 ± 0.007 quoted in Ref. The quantity z^ is a microscopic ('extrapolation') 
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length depending on Ji/J, which ensures that mi = m(0) is finite. It drops out in the 
scahng regime z ^ Ze- In accordance with our expectations we found it to have the same 
value for T > Tc and T < Tc. For large distances z, which are still sufficiently small so 
that finite size effects are negligible, the profiles ■m±{z) decay exponentially on the scale of 
^, where the value of ^ is in conformity with the results of Ref. |T^. At even larger values 
of z finite size effects are clearly visible. 

In order to check whether scaling holds we analyzed the data for m{z) in the first twenty 
layers closest to the surface, using the scaling ansatz 

m{z) ~ M^\TfP±{^) (19) 

for reduced temperatures r in the ranges given in Figs. 1 and 2. As is borne out by these 
figures, the data collapsing on a single curve P± works reasonably well both for r > 
and r < 0, within the limitations caused by statistical error and finite size effects. The 
resulting Monte Carlo profiles -P±(C) agree qualitatively quite well with our extrapolated 
e expansion results. Quantitatively, the agreement is not very satisfactory. To appreciate 
these results, one should note, however, that our extrapolated RG results are free of any 
adjustable parameter. Likewise, no adjustable parameter remains in the Monte Carlo results 
for P±, since the scales and amplitudes are fixed by the required large-distance forms ([12D 
and (plsl) . (One probably could get improved agreement at intermediate values of ( by 
exploiting the freedom to make O(e^) errors at order e, but in contrast to our procedure, 
this would be unsystematic.) 

The results for P± may be used to estimate various universal ratios. The ratio 

Rma ^ / d(P^{C)/ / dC [P_(C) - 1] , (20) 
JO Jo 

which is proportional to the ratio of the amplitudes of the excess order parameter above 



and below Tc, has been investigated experimentally by Law and Smith |T6| recently. Integra- 
tion of our Monte Carlo results yields Rma — 1-14 in good agreement with their experimental 
estimate 1.18 ± 0.13. Integration of our extrapolated RG results yields a somewhat higher 
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value ~ 1.3. Of considerable interest are also the universal ratios c±. From our Monte Carlo 



results we find the estimates c+ = 0.866±0.07 and c_ = 1.22±0.08. Floter and Dietrich ||T2 
estimated c+ theoretically by means of Fade approximants into which both the e expansion 
results of Ref. ||^ and the exact results for d = 2 were incorporated. They also extracted 
values of c+ from a number of experiments. These experimental estimates were found to 
scatter around their theoretical estimate c+ ~ 0.95, which is in accordance with our Monte 
Carlo estimate. Finally, we mention that the universal ratio P'^i/cj^ was considered first 
by Liu and Fisher [Q. From their analysis of critical adsorption data they obtained the 
estimate ^ 0.85. Our Monte Carlo data yield a considerably larger value of ^ 1.85. This 
discrepancy might be due to the limited accuracy of the experimental data available at that 
time. A problem could also be that the approximate forms of -P+(C) used in the analysis of 
the experimental data were not accurate enough. (As discussed by Liu and Fisher 0] the 
estimate is rather sensitive to the chosen form of P+.) 

4. CONCLUDING REMARKS 



In summary, we presented results for the universal scaled order parameter profiles P±{C) 
at critical adsorption and the extraordinary transition. These were obtained by two distinct 
methods, namely, by RG-improved perturbation theory in 4 — e dimensions and by Monte 
Carlo simulations. To obtain meaningful extrapolations of the RG results to the physically 
interesting dimension = 3, a novel renormalization scheme was used that guarantees the 
correct exponentiation of the leading short- distance singularity. Both the latter method and 
the Monte Carlo simulations yield profiles for (i = 3, which exhibit the expected smooth 
and monotonous crossover from the power-law behavior (|12D at small C to the large-distance 
form ([T^). The results of both methods are in fair agreement with each other. The knowl- 
edge of the profiles can be utilized to estimate various universal amplitude ratios, which 
can also be determined by experiments. The agreement of our theoretical estimates with 
the few available experimental estimates ITOIT^I of the universal numbers c+ and Rma is 
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encouraging. Yet, much more detailed experimental checks of the theoretical predictions are 
clearly needed. 
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FIGURES 



Fig. 1: Scaling function P+{C). The dashed curve shows the mean-field approximation. 
The full curve is the extrapolation to d 3 of the one-loop RG result obtained with the help 
of our specially adapted renormalization scheme. The symbols show Monte Carlo results for 
different reduced temperatures. 

Fig. 2: Scaling function -P-(C)- The dashed and full curves are the r < analogs 
of the curves in Fig. 1. The symbols represent Monte Carlo results for different reduced 
temperatures. 
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Fig. 1: Scaling function -P+(C)- The dashed curve shows the mean-field approximation. 
The full curve is the extrapolation to d = 3 of the one-loop RG result obtained with the help 
of our specially adapted renormalization scheme. The symbols show Monte Carlo results for 
different reduced temperatures. 
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Fig. 2: Scaling function -P-(C)- The dashed and full curves are the r < analogs 
of the curves in Fig. 1. The symbols represent Monte Carlo results for different reduced 
temperatures. 
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